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We propose a modification in the Sznajd sociophysics model defined on the square lattice. For this 
purpose, we consider reputationa mechanism limiting the agents persuasive power. The reputation 
is introduced as a time-dependent score, which can be positive or negative. This mechanism avoids 
dictatorship (full consensus, all spins parallel) for a wide range of model parameters. We consider two 
different situations: case 1, in which the agents reputation increases for each persuaded neighbor, and 
case 2, in which the agents reputation increases for each persuasion and decreases when a neighbor 
keeps his opinion. Our results show that the introduction of reputation avoids full consensus even 
for initial densities of up spins greater than 1/2. The relaxation times follow a log- normal- like 
distribution in both cases, but they are larger in case 2 due to the competition among reputations. 
In addition, we show that the usual phase transition occurs and depends on the initial concentration 
d of individuals with the same opinion, but the critical points d c in the two cases are different. 

Keywords: Dynamics of Social Systems, Phase Transitions, Cellular Automata. 



I. INTRODUCTION 

Ising-type models have been reviewed and used by 
physicists in many different areas, such as sociology, pol- 
itics, marketing, and finance [lHj]. I n 2000, an agent- 
based model proposed by Sznajd- Weron et al. [f_] was 
successfully applied to the dynamics of a social system. 
In particular, the model reproduced certain properties 
observed in a real community. At the focus of the Sz- 
najd model (SM) is the emergence of social collective 
(macroscopic) behavior due to the interactions among 
individuals, which constitute the microscopic level of a 
social system. 

This model has been extensively studied since the in- 
troduction of the original one-dimensional model in 2000. 
Modifications were proposed in many works, like the con- 
sideration of different types of lattices such as square [_| , 
triangular [6j and cubic [7[, the increase of the range of 
the interaction [8fl and the number of variable's states 
[9l-[ll| and the possibility of diffusion of agents [ill EH ■ 

The original SM consists of a chain of sites with pe- 
riodic boundary conditions where each site (individual 
opinion) could have two possible states (opinions) repre- 
sented in the model by Ising spins ("y es " or "no"). A 
pair of parallel spins on sites i and i + 1 forces its two 
neighbors, i — 1 and i + 2, to have the same orientation 
(opinion), while for an antiparallel pair (i, i + 1), the left 
neighbor (i — 1) takes the opinion of the spin i + 1 and the 
right neighbor (i + 2) takes the opinion of the spin i. In 
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this first formulation of the SM two types of steady states 
are always reached: complete consensus (ferromagnetic 
state) or stalemate (anti- ferromagnetic state), in which 
every site has an opinion that is different from the opin- 
ion of its neighbors. However, the transient displays a 
interesting behavior, as pointed by Stauffer et al. 
Defining the model on a square lattice, the authors in 
considered not a pair of neighbors, but a 2 x 2 plaque- 
tte with four neighbors. Considering that each plaquette 
with all spins parallel can convince all their eight neigh- 
bors (we will call this StaufFer's rule), a phase transition 
was found for an initial density of up spins d = 1/2. 

It is more realistic to associate a probability of persua- 
sion to each site. The SM is robust with respect to this 
choice: if one convinces the neighbors only with some 
probability p, and leaves them unchanged with probabil- 
ity 1 — p, still a consensus is reached after a long time 
Q]. Models that consider many different opinions (us- 
ing Potts' spins, for example) or defined on small- world 
networks were studied in order to represent better ap- 
proximations of real communities' behavior (see Q and 
references therein). In another work, in order to avoid 
full consensus in the system and makes the model more 
realistic, Schneider introduced opportunists and persons 
in opposition, that are unconvinced by their neighbors 

a 

In the real world, however, the dynamics of social rela- 
tionships is more complex. Even when such more struc- 
tured topologies as small-world networks are adopted to 
bring the SM closer to reality, a large number of details 
is often neglected. In order to advance toward realism, 
we recently considered a reputation mechanism [14j . We 
believe that the reputation of agents who hold the same 
opinion is an important factor in opinion propagation 
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across the community. In other words, it is realistic to 
believe that the individuals will change their opinions un- 
der the influence of highly respected persons. The repu- 
tation limits the agents power of persuasion, and we can 
expect the model in 1J] to be more realistic than the 
standard one In fact, we showed that simple micro- 
scopic rules are sufficient to generate a democracy-like 
state, ferromagnetically ordered with only partial polar- 
ization [i4| . 

In this work, we revise and extend our previous re- 
sults by allowing reputations to increase and decrease, 
depending on whether the agents are or are not per- 
suaded. This generalization is based on the behavior of 
real social networks: Certain persons tend to be skeptical 
if the persuaders have low reputation, in which case their 
best strategy is to keep their opinions. In this sense, in- 
cluding reputation makes the SM more realistic. To be 
thorough, we will consider two different protocols. In the 
first case, the agents' reputations increase for each per- 
suaded neighbor, whereas in the second case, the agents' 
reputations rise in case of persuasion and decrease when- 
ever the agents fail to convince their neighbors. 

This paper is organized as follows: in Section |TT] we 
present the model and define their microscopic rules. The 
numerical results as well as the finite-size scaling analysis 
are discussed in Section IIIII Finally, in Section IIVI we 
summarize our conclusions. 



II. MODEL 

We have considered our model defined on a square 
lattice with L x L agents and periodic boundary condi- 
tions. Similar to Stauffer's rule (rule la of |4|), we choose 
at random a 2 x 2 plaquette of four neighbors and if 
all central spins are parallel, the neighbors may change 
their opinions. The difference in our model is that the 
neighbor's spins will be flipped depending on the pla- 
quette reputation. An integer number (R) labels each 
player and represents its reputation across the commu- 
nity, in analogy to the Naming game model considered 
by Brigatti [15] . The reputation is introduced as a score 
for each player and is time dependent. The agents start 
with a random distribution of R values, and during the 
time evolution, the reputation of each agent changes ac- 
cording to its capacity of persuasion. We will consider 
in this work that the initial values of the agents' repu- 
tation follow a gaussian distribution centered at with 
standard deviation a. 

We have considered two different situations in this 
work: in the first case, the reputations increase following 
the model's rules, whereas in the second case the 
reputations may increase and decrease. One time step 
in our model is defined by the following microscopic rules: 



1. We randomly choose a 2 x 2 plaquette of four 
neighbors; 

2. If not all four center spins are parallel, leaves its 
eight neighbors unchanged. 

3. On the other hand, if the four center spins are fully 
polarized, we calculate the average reputation R of 
the plaquette, 



R 



zX>. 
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where each term Ri represents the reputation of 
one of plaquettes' agent. 

4. We compare the reputations of each of the eight 
neighbors of the plaquette with the average repu- 
tation. If the reputation of a neighbor is less than 
the average, this neighbor follow the plaquette ori- 
entation. On the other hand, if the neighbor's rep- 
utation exceeds R, no action is taken. 

5. For each persuasion, the reputation of the plaque- 
tte agents is incremented by 1, so that the average 
plaquette reputation is increased by 1. 

Case 2 

In this case, steps 1-4 are as described above. Step 
5, by contrast, is changed to the following rule: 

• For each persuasion, the reputation of the plaquette 
agents is incremented by 1. On the other hand, for 
each failure, the reputations within the plaquette 
are decremented by 1. 

Thus, even in the case of fully polarized plaquettes, 
different numbers of agents may be convinced, namely 
8, 7, 6, . . . , 1, or 0. As pointed by Stauffer in 
[2|, we can imagine that each agent in the Sznajd model 
carries an opinion that can either be up (e.g., Republican) 
or down (e.g., Democrat), which represents one of two 
possible opinions on any question. The objective of the 
agents in the game is to convince their neighbors. One 
can expect that, if a certain group of agents convince 
many others, their persuasive power grows. On the other 
hand, the persuasive powers may drop if the agents fail to 
convince other individuals. The inclusion of reputation 
in our model captures this feature of the real world. 



III. NUMERICAL RESULTS 



A. Case 1: Emergence of consensus 



Case 1 



In the simulations, we considered er = 5. Following 
the previous works on the SM, we can start studying the 
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FIG. 1. Time evolution of the magnetization (case 1) for 
L = 53, initial densities of up spins d = 0.4 and d = 0.6 
and different samples (a). We can see that the steady states 
show situations where the total consensus is not obtained, in 
opposition of the standard Sznajd model defined on the square 
lattice [4|. In figure (b) we show the results for d — 0.1 and 
d = 0.9. In these cases the system reaches consensus in all 
samples. 



time evolution of the magnetization per site, 

N 



7)1 



(i) 



where N = L 2 is the total number of agents and Sj = ±1. 
In the standard SM defined on the square lattice Q, the 
application of the Stauffer's rule, where a 2 x 2 plaquette 
with all spins parellel convince its eigth neighbors, with 
initial density of up spins d = 1/2 leads the system to 
the fixed points with all up or all down spins with equal 
probability. For d < 1/2 (> 1/2) the system goes to a 
ferromagnetic state with all spins down (up) in all sam- 
ples, which characterizes a phase transition at d = 1/2 
in the limit of large L. As pointed by the authors in [J], 
fixed points with all spins parallel describe the published 
opinion in a dictatorship, which is not a commom situ- 
ation nowadays. However, ferromagnetism with not all 
spins parallel corresponds to a democracy, which is very 
commom in our world. We show in Fig. [T]the behavior 
of the magnetization as a function of the simulation time 
in our model, for case 1. In Fig. [T](a), we show a value of 
d>l/2(<l/2), and one can see that the total consensus 
with all spins up (down) will not be achieved in any sam- 
ple. On the other hand, in Fig. [1] (b) we show situations 



FIG. 2. Log-log plot of the histogram of relaxation times (case 
1) for L = 53 and d — 0.8, obtained from 10 4 samples, with 
agents' initial reputations following a gaussian distribution 
with different standard deviations a (a). The distribution is 
compatible with a log-normal one for all values of a, which 
corresponds to the observed parabola in the log-log plot. It is 
also shown the average relaxation time r, over 10 4 samples, 
versus latice size L in the log- log scale (b). The straight line 
has slope 5/2. The result is robust with respect to the choice 
of different a values. 



where the consensus is obtained with all up (for d = 0.9) 
and all down spins (for d = 0.1). These results indi- 
cate that (i) a democracy-like situation is possible in the 
model without the consideration of a mixing of different 
rules [J], or some kind of special agents, like contrarians 
and opportunists and (ii) if a phase transition also 
occurs in our case, the transition point will be located 
somewhere at d > 1/2. 

We have also studied the relaxation times of the model, 
i.e., the time needed to find all the agents at the end hav- 
ing the same opinion. The distribution of the number of 
sweeps through the lattice, averaged over 10 4 samples, 
needed to reach the fixed point is shown in Fig. [2] (a). 
We can see that the relaxation time distribution is com- 
patible with a log-normal one for all values of the stan- 
dard deviation a, which corresponds to a parabola in the 
log-log plot of Fig. [2] (a). The same behavior was ob- 
served in other studies of the SM d [H \]]% . In Fig. H 
(b) we show the average relaxation time r [also over 10 4 
samples, considering the relaxation times of Fig. [2] (a)] 
versus latice size L in the log-log scale. We can verify a 
power-law relation between these quantities in the form 
r ~ L 5 / 2 , for large L and all values of the standard de- 
viation, which indicates that this result is robust with 
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FIG. 3. Fraction / of samples (case 1) which show all spins 
up when the initial density of up spins d is varied in the range 
0.4 < d < 1.0, for some lattice sizes L (a). The total number 
of samples are 1000 (for L = 31 and 53), 500 (for L — 73 and 
101) and 200 (for L = 121). It is also shown the corresponding 
scaling plot of / (b). The best collapse of data was obtained 
for a = 0.03, b = 0.47 and d c = 0.88. 



respect to the choice of different a values. Power-law re- 
lations between r and L were also found in a previous 
work on the SM [l]|. 

We can now analyze the phase transition of the model. 
For this purpose, we have simulated the system for dif- 
ferent lattice sizes L and we have measured the fraction 
of samples which show all spins up when the initial den- 
sity of up spins d is varied in the range 0.4 < d < 1.0. 
In other words, this quantity / give us the probability 
that the population reaches consensus, for a given value 
of d. We have considered 1000 samples for L = 31 and 
53, 500 samples for L = 73 and 101 and 200 samples 
for L = 121. The results are shown in Fig. [3] (a). One 
can see that the transition point is located somewhere 
in the region d > 1/2, as above discussed. In order to 
locate the critical point, we performed a finite-size scal- 
ing (FSS) analysis, based on the standard FSS equations 

mm 



f(d,L) = L~ a f((d-d c )L b ) , 
d c (L) = d c + c L- b , 



(2) 
(3) 



where c is a constant and / is a scaling function. The 
result is shown in Fig. [3] (b) , and we have found that 



d c = 0.88 ±0.01 



(4) 



FIG. 4. Fraction / of samples (case 1) which show all spins up 
when the initial density d is varied in the range 0.4 < d < 1.0, 
for L — 53, 1000 samples and some different values of a. This 
result show that the increase of a do not change the behavior 
of /. 



in the limit of large L. In addition, we have obtained 
a = 0.030 ± 0.005 and b = 0.47 ± 0.02. The critical point 
occurs at d > 1/2, different of the SM without reputa- 
tion defined on the square lattice. This fact may be eas- 
ily understood: at each time step, the randomly choosen 
2x2 plaquette may convince 8, 7, 6, 1 or neighbors, 
even if the plaquettes' spins are parallel. In the standard 
model, if the plaquettes spins' orientations are the same, 
8 neighbors are convinced immediately, thus it is neces- 
sary a smaller initial density of up spins to the system 
reaches the fixed point with all spins up. Thus, the usual 
phase transition of the SM also occurs in our model, in 
case 1, and this transition is robust with respect to the 
choice of different values of a (see Fig. [4j . 



B. Case 2: Competition among reputations 

As discussed in section [Til m this second case the 
agent's reputations may increase and decrease, which de- 
fines a competition of reputations in the game. The evo- 
lution of the magnetization per site is shown in Fig. [S] 
In the case of intermediary densities d the system reaches 
steady states with m < 1, i.e., we have democracy-like 
situations. However, due to the competition of reputa- 
tions, that increase and decrease depending on the av- 
erage reputation of the plaquettes during the time evo- 
lution, the system reaches steady states with different 
magnetizations. Another consequence of the competi- 
tion appears in the case of large and very small initial 
densities d: even for the cases d = 0.9 and d = 0.1 the 
system reaches consensus only in some realizations of the 
dynamics. This fact can be observed in the inset of Fig. 
[S](b): the dotted line is m = 1, and we observe that just 
one of the three realizations reaches consensus. Thus, 
for case 2, the consensus is very hard to be obtained. 
Nonetheless, the emergence of democratic steady states 
is favoured in this second case, in comparison with case 
1, which makes the model more realistic in this sense. 
We have also studied the relaxation times for case 2. 
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FIG. 5. Time evolution of the magnetization (case 2) for 
L = 53, initial densities of up spins d = 0.4 and d = 0.6 and 
different samples (a). We can see differences between these 
steady states and those of case 1, but we also have democracy- 
like situations. In figure (b) we show the results for d — 0.1 
and d = 0.9. Observe in the inset that even for large values 
of densities like d = 0.9 the system reaches consensus only in 
some samples (analogously for d = 0.1). The dotted line in 
the inset is m = 1 (full consensus). 



FIG. 6. Log-log plot of the histogram of relaxation times (case 
2) for L = 53 and d = 0.99, obtained from 10 4 samples, with 
agents' initial reputations following a gaussian distribution 
with different standard deviations a (a). The distribution is 
compatible with a log-normal one for all values of a, which 
corresponds to the observed parabola in the log-log plot. It is 
also shown the average relaxation time r, over 10 4 samples, 
versus latice size L in the log-log scale (b). The power-law 
behavior for large L is r ~ L 3 , for all values of a. 



The distribution of the number of sweeps through the 
lattice, averaged over 10 4 samples, needed to reach the 
fixed point is shown in Fig. [6] (a) for some values of 
the standard deviation a. We can see that, as in case 
1, the relaxation time distribution is compatible with a 
log-normal one for all values of a. However, due to com- 
petition of reputations, the relaxation times of case 2 are 
greater than the corresponding relaxation times of case 
1. In Fig. [6] (b) we show the average relaxation time r 
[also over 10 samples, considering the relaxation times 
of Fig. [6] (a)] versus latice size L in the log-log scale. 
In this case, we verify the power-law behavior t ~ L 3 
for large L and all values of the standard deviation. In 
other words, the competition of reputations increases the 
relaxation times of the system, as above discussed, and 
this effect becomes stronger when we increase the num- 
ber of agents of the system (or the lattice size L), which 
implies in a power-law exponent for r greater than the 
exponent for the case 1. 

Following the approach of the last subsection (case 1), 
we have simulated the system for different lattice sizes 
L and we have measured the fraction of samples which 
show all spins up when the initial density of up spins d 
is varied in the range 0.5 < d < 1.0. We have considered 



the same number of samples of the last subsection, and 
the results are shown in Fig. [7] (a). One can see that the 
transition point is located somewhere in the region d > 
0.88, i.e., the critical density in case 2 is greater than in 
case 1, as expected due to the competition of reputations. 
We determined the critical point for this case using the 
above Eqs. © and ©• The best collapse of data is 
shown in Fig. [7J (b), obtained with d c = 1.00 ± 0.01, 
a = 0.00 ± 0.01 and b = 1.37 ± 0.02. In other words, the 
case 2 presents a different critical density and different 
critical exponents, in comparison with case 1. However, 
the usual phase transition of the SM also occurs in case 
2, and this transition is robust with respect to the choice 
of different values of a (see Fig. [8]). 

Observe that, in order to minimize the finite-size ef- 
fects, we excluded the smaller size L — 31 for the FSS 
process [see the inset of Fig. [7J (b)]. In fact, we can ob- 
serve in Fig. [7J (a) that, for L = 31, the curve of the 
quantity / presents a inflection point, which not appears 
in the other sizes. Thus, this inflection point in the curve 
for L = 31 is a pronounced finite-size effect of the model 
considered in case 2. 
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FIG. 7. Fraction / of samples (case 2) which show all spins 
up when the initial density of up spins d is varied in the range 
0.5 < d < 1.0, for some lattice sizes L (a). The total number 
of samples are 1000 (for L = 31 and 53), 500 (for L — 73 and 
101) and 200 (for L = 121). It is also shown the corresponding 
scaling plot of / (b). The best collapse of data was obtained 
for a — 0.0, 6 = 1.37 and d c — 1.0. To minimize the finite- 
size effects, we have excluded the smaller size L = 31 of the 
collapse (see the inset). 
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FIG. 8. Fraction / of samples (case 2) which show all spins up 
when the initial density d is varied in the range 0.5 < d < 1.0, 
for L = 53, 1000 samples and some different values of a. This 
result show that the increase of a do not change the behavior 
of/. 



introduced as a score for each player and is time depen- 
dent, varying according to the model's rules. The agents 
start with a random distribution of reputation values, 
and during the time evolution, the reputation of each 
agent changes according to its capacity of persuasion. 
We have considered in this work that the initial values 
of the agents' reputation follow a gaussian distribution 
centered at with standard deviation a. In addition, we 
have studied separately two different situations: (i) a case 
where the reputations increase due to each persuaded in- 
dividual (case 1), and (ii) a case where the reputations 
increase for persuasion and decrease if a group of agents 
fail to convince one of its neighbors (case 2). 

In the first case, we observed a log-normal-like distri- 
bution of the relaxation times, i.e., the time needed to 
find all the agents at the end having the same opinion. 
In addition, the average relaxation times grow with the 
linear dimension of the lattice in the form r ~ L 5 / 2 . The 
system undergoes the usual phase transition, that was 
identified by measurements of the fraction / of samples 
which show all spins up when the initial density of up 
spins d is varied. In other words, this quantity / give 
us the probability that the population reaches consen- 
sus, for a given value of d. We localized the transition 
point by means of a finite-size scaling analysis, and we 
found d c — 0.88. This critical density is greater than 
1/2, the value found by Stauffer et at [J] in the standard 
formulation of the Sznajd model. This fact may be eas- 
ily understood: at each time step, the randomly choosen 
2x2 plaquette may convince 8, 7, 6, 1 or neighbors, 
even if the plaquettes' spins are parallel. In the standard 
case, if the plaquettes spins' orientations are the same, 8 
neighbors are convinced immediately, thus it is necessary 
a smaller initial density of up spins to the system reaches 
the fixed point with all spins up. The simulations in- 
dicate that the observed phase transition is robust with 
respect to the choice of different values of a. 

In the second case, the steady states with m < 1 are 
favoured due to the competition of reputations, and even 
for large densities d the system reaches consensus only in 
some samples. We also found that the relaxation times 
are log-normally distributed, but they are greater than 
the relaxation times of case 1. The average relaxation 
times are greater than the corresponding values found in 
the first case, and they also depend on the linear lattice 
size in a power-law form, r ~ L 3 , but with a greater 
exponent in comparison with case 1. The usual phase 
transition also occurs in case 2, but the critical density 
was found to be d c = 1.0. In addition, the second sit- 
uation presents strong finite-size effects. The observed 
differences between the two cases are due to the compe- 
tition of reputations that occurs in case 2. 



IV. CONCLUSIONS 



In this work we have studied a modified version of 
the Sznajd sociophysics model. In particular we have 
considered reputation, a mechanism that limits the ca- 
pacity of persuasion of the agents. The reputation is 
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